o 
o 

(N 

X> 

IX, 



Multiple-pulse coherence enhancement of solid state spin qubits 
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We describe how the spin coherence time of a localized electron spin in solids, i.e. a solid state spin 
qubit, can be prolonged by applying designed electron spin resonance pulse sequences. In particular, 
the spin echo decay due to the spectral diffusion of the electron spin resonance frequency induced 
by the non-Markovian temporal fluctuations of the nuclear spin flip-flop dynamics can be strongly 
suppressed using multiple-pulse sequences akin to the Carr-Purcell-Meiboom-Gill pulse sequence 
in nuclear magnetic resonance. Spin coherence time can be enhanced by factors of 4-10 in GaAs 
quantum dot and Si:P quantum computer architectures using composite sequences with an even 
number of pulses. 

PACS numbers: 03.67.-a; 76.60.Lz: 03.65.Yz: 76.30.-v; 03.67.Lx 
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Understanding spin coherence in solids is one of the 
oldest problems in condensed matter physics [l|, |2| , go- 
ing back to the seminal pioneering work of Hahn [3] , Fe- 
her [4J, and others. Recently, there has been a resurgence 
of widespread interest in this problem in the context 
of quantum computation using spin qubits in semicon- 
ductors. In particular, a necessary condition for quan- 
tum computation is long qubit coherence time so that 
quantum error correction can be meaningfully imple- 
mented. For various proposed electron spin qubits in 
semiconductor nanostructures, coherence times on the 
order of 0.1 ms (or longer) are required as a necessary 
condition; assuming typical gating times of 100 ps, this 
leads to a Q-factor of ~ 0.1 ms/100 ps = 10 6 which 
satisfies the current quantum error correction constraint 
(~ 10~ 4 — 10~ 6 ). In fact, the only reason for the current 
experimental and theoretical interest in spin-qubit based 
scalable solid state quantum computer architectures (as 
opposed to charge qubits, which are easier to manipulate 
and read out) is the expected long spin coherence times 
G^ M s ) compared with charge (i.e. orbital) coherence 
times (< ns) in solids, making it inrpossible (possible) 
for quantum error correction schemes to work for solid 
state charge (spin) qubits. Understanding all aspects of 
electron spin decoherence in solid state nuclear spin en- 
vironment is important so that effective strategies can 
be developed to restore and/or enhance coherence by, for 
example, efficient pulse engineering. 

It is, therefore, highly desirable to enhance the co- 
herence time of spin qubits in semiconductor struc- 
tures, particularly since semiconductor-based spin quan- 
tum conrputer architectures have considerable advan- 
tages in terms of scalability and fabrication. It has been 
known for a long time that various refocusing techniques 
using electromagnetic pulses (e.g. Hahn spin echo) sub- 
stantially enhance spin coherence compared with the free 
induction decay. In this Letter we describe techniques 
that substantially enhance coherence of solid state spin 
qubits by utilizing various composite sequences going be- 



yond the simple Hahn spin echo. Demonstrating suffi- 
cient spin coherence in such a way is a necessary step 
in the development of spin-based fault tolerant quantum 
computing. This Letter has direct implications for quan- 
tum spin memory in semiconductors and offers potential 
for coherent quantum logic gate design. 

A particular impetus for our theoretical study comes 
from the beautiful recent experimental work by Petta et 
al. [5| who measured Hahn spin echo decay of electron 
spin qubits in GaAs coupled quantum dot architectures, 
finding a spin coherence time Ti ~ 1 /lis (consistent with 
our spectral diffusion theory of spin decoherence in GaAs 
quantunr dot spin qubits [6j), substantially enhancing (by 
a factor of 100) the inhomogeneous spin dephasing time 
of T 2 * ~ 10 ns in the same system. In the current Letter, 
we investigate the potential for further coherence time 
enhancement by utilizing a more complex pulse sequence 
originally used by Carr and Purcell a long time ago [7| 
in the context of magnetic resonance studies. Further 
coherence enhancement in GaAs quantum dots has re- 
cently been achieved using the multiple pulse sequence 
suggested in this Letter, and the experimental results 
are in good agreement with our predictions 8]. Earlier 
experinrental work on the ensemble of P donor electronic 
states in Si, of great interest to Si quantum computer 
architectures, reported very long (> ms) spin echo Ti 
coherence times [9J, [lfj] , which can be further enhanced 
in Si through isotopic purification. We find that com- 
posite sequences beyond the simple Hahn echo sequence 
could lead to substantial spin coherence enhancement in 
Si spin qubits also. 

We consider a single localized electron spin in a semi- 
conductor (e.g., GaAs quantum dot, Si:P) interacting 
with the surrounding nuclear spin bath. The electron's 
spin polarization is preserved by applying a magnetic 
field due to the very large (a factor of ~ 2000) differ- 
ence between the electron and nuclear Zeeman energy 
splittings. Virtual electron spin-flip transitions allow 
hyperfine-mediated interactions between nuclei, but such 



processes cause only a small visibility decay of refocused 
echoes that is sufficiently small for fault tolerant quan- 
tum computing at fields above B ~ 1 T [ll| in GaAs 
quantum dots for which the effect is particularly strong. 
Echo modulation due to anisotropic hypcrfine coupling 
is also sufficiently suppressed at fields above B ~ 9 T 
[12J ] in Si:P for which the effect is particularly strong. 
In such a situation (i.e. an applied magnetic field and 
low < 100 mK operational temperatures), it is now well 
accepted that the decoherence of a solid state electron 
spin qubit is dominated by the spectral diffusion process, 
where the flip-flop dynamics of interacting nuclear spins 
creates a temporally random non-Markovian magnetic 
field at the electron spin location leading to decoherence. 
We have earlier theoretically studied spectral diffusion ef- 
fects on the Hahn spin echo introducing a quantum clus- 
ter expansion technique [y, [l3| that provides a formally 
exact treatment of the non-Markovian quantum nuclear 
dynamics in the large applied magnetic field limit. In this 
Letter we apply this technique to multi-pulse echo se- 
quences in order to investigate the possible enhancement 
of spin coherence. We find that carefully designed multi- 
pulse sequences could considerably enhance solid state 
spin qubit coherence with the spin decoherence time T^ 
[14| increasing strongly with the number of pulses. Our 
work has obvious implications for the design and opera- 
tion of spin quantum computer architectures in semicon- 
ductor nanostructures. In particular, we find that our 
non-Markovian treatment leads to low order symmetry 
related cancellations when an even number of pulses are 
applied. We also predict that the logarithm of the echo, 
as a function of inter-pulse time (i.e., r), scales as the 
square of the number of applied pulses in stark contrast 
to linear scaling of the stochastic theory [15| . 

Spins of interest (e.g., spins representing qubits) can 
be refocused, partially reversing the effects of their local 
magnetic fields, by applying to them sequences of rotat- 
ing pulses (e.g., via resonance). Our analysis treats these 
pulses as ideal, making the approximating assumption 
that they perform exact rotations of some desired spin 
in an infinitesimal time without affecting the rest of the 
system. We consider the application of the Carr-Purcell- 
Meiboom-Gill (CPMG) [7(] pulse sequence (equivalent, 
for our purposes with ideal pulses, to the Carr-Purcell 
sequence) to a solid state electron spin qubit. We repre- 
sent this sequence as (r — ► it — > r) n by which we mean 
that the system evolves freely for a time t, we apply a 
7r-rotation pulse perpendicular to the applied magnetic 
field, evolve the system freely for time r again, and re- 
peat this process n times. The n = 1 case is equivalent 
to the Hahn echo. Except for comparison with the Hahn 
echo, we treat the n = 1v case of an even number of 
pulses. We find analytically (within relevant approxima- 
tions) that even pulses are enhanced. A very general 
cluster expansion technique was introduced in Ref. [13l |, 
where it was successfully applied to Hahn echoes of Si:P 



and was recently verified in an independently developed 
theory [l6(. We apply this technique now to the CPMG 
sequence. This analysis begins by writing the exact ex- 
pression for the echo of the (t — > it 



t) 2l/ sequence: 
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where the first summation is re- 
stricted to pairs of like nuclei (so that Zeeman energy 
is preserved when they flip-flop). The Tr operation in 
Eq. {T]) traces over the states of the nuclear bath and 
M is the number of such states. The {A n } are hyper- 
fine coupling constants between the spins of the electron 
and nucleus n; typically max (A n ) ~ 10 6 s _1 (with h = 1 
units) . Nucleus n has spin operators denoted with I n and 
its gyromagnetic constant is j n . The {b nm } are dipolar 
coupling constants between nuclei in the bath; typically, 
max(6 nm ) ~ 10 2 s" 1 so that b nm <C A n which is impor- 
tant for the convergence of our cluster expansion. Details 
are given in Ref. 6] . 

Our cluster expansion is an approximation that is 
based upon the fact that we may exactly decompose 
Eq. ^ into a sum of all possible products of contribu- 
tions from disjoint sets of nuclei: 
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The full proof and necessary conditions for such a de- 
composition is provided in Ref. 6]. Each Ci denotes a 
set of nuclei and v' c .(t) is called the "contribution" from 
this set. A contribution has the property that it can only 
be significant when interactions between nuclei in the set 
are significant and no part is isolated from the rest. We 
consider only local dipolar interactions in the current Let- 
ter, and thus contributions only arise when the nuclei in 
the set are spatially clustered together; hence we refer to 
these sets as clusters. Before we define the "cluster con- 
tribution," v' c (t), we first define Ue(r) as the solution to 
Eq. |T]) when only considering nuclei contained in some 
set (or cluster) C. A cluster contribution may then be 
recursively defined by and computed using 



v 'c( T ) = v c(t) 
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{Ci} disjoint. 
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subtracting from vq (t) the sum of all products of contri- 
butions from disjoint sets of clusters contained in C. 

The decomposition of Eq. ((2]) is useful when making 
an approximation that assumes cluster contributions de- 
crease with an increase in cluster size. We can then define 
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FIG. 1: (color online). Numerical results for echo envelopes 
of the Hahn and various CPMG pulse sequences (labelled by 
number of pulses, 2v) for a quantum dot in GaAs with a 50 nm 
Fock-Darwin radius, 10 nm quantum well thickness, and ap- 
plied magnetic field and thickness direction both along [001] 
in the lattice. Solid lines show convergent cluster expansion 
(exact) results while dotted lines show lowest order perturba- 
tive results (within the cluster expansion framework). 




FIG. 2: (color online). Numerical results for echo envelopes 
of the Hahn and various CPMG pulse sequences (labelled by 
number of pulses, 2v) for the Si:P donor electron in natural 
Si with the applied magnetic field along the [001] lattice di- 
rection. Solid lines show convergent cluster expansion (exact) 
results while dotted lines show lowest order perturbative re- 
sults (within the cluster expansion framework). The ordinate 
range is the first 10% (0.9 to 1) of echo decay where the cluster 
expansion is convergent for all pulse sequences shown. 
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Because of the difficulty of iterating through all possi- 
ble disjoint sets of clusters (below a given size), we make 
a further approximation by relaxing the constraint that 
they be disjoint sets. In this approximation, the loga- 
rithm of the fc th order result simply becomes the sum of 
all cluster contributions up to size k: 
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This approximation may be tested (or corrections made) 
as described in Ref. [g|. We wish to note the simplicity 
of implementing our approach using Eqs. ^ and @ as 
compared to a very recent diagrammatic approach [17| ] . 
We show numerical results obtained by this method 
(up to visible convergence) in Figs. [T] and [2] for the coher- 
ence of a quantum dot electron in GaAs and for the coher- 
ence of an electron bound to a donor in Si:P respectively, 
comparing Hahn and CPMG echo envelopes Hahn echo 
envelopes. The effect of the coherence enhancement for 
even echoes is noticeable in both systems by comparing 
results of two-pulse CPMG versus the Hahn echo. These 
figures also show that coherence is indeed prolonged with 
an increase in the number of pulses. We note that these 
are plotted with respect to the total decay time, t = Avt 
(for the Hahn echo, t — 2r). 



There are two perturbation theories that help to ex- 
plain cluster expansion convergence via diminishing clus- 
ter contributions with increasing cluster size: the "inter- 
bath perturbation" (previously called dipolar perturba- 
tion [6() that treats the coupling between nuclei in the 
bath as a small perturbation relative to the (potentially) 
strong coupling of the nuclei to the electron, and the 
time (or r) expansion. Using A to generically denote the 
(small) perturbation parameter of either theory, it was 
shown [6] that a contribution from a cluster of size k is 
0(X k ). The inter-bath perturbation is most applicable 
to nuclei that are coupled strongly to the electron (e.g., 
those closer to the center), and the time expansion is 
most applicable to nuclei with a weak coupling to the 
electron (e.g., those further from the nuclei where lower 
energies result in slower evolution). The two theories 
work in tandem to produce overall convergence of the 
cluster expansion for a variety of physical systems. 

By using the lowest order results of cluster contribu- 
tions [Eq. ([3])] with respect to either perturbation theory 
in the computation of the cluster expansion [Eq. (JS}] , we 
can see the important role played by these perturbation 
theories and gain useful insights. The lowest order result 
with respect to inter-bath perturbation explains the en- 
hancement of even CPMG echoes relative to odd echoes 
as well as the scaling behavior with the number of pulses. 
Using A ~ b nm /\A n — A m \ as the perturbation parameter, 
the low order solution to Eq. (Q]) applied to a given clus- 
ter, vc(t), behaves as ve(r) = 1 — 0(^ 2 A 4 ). For a cluster 
contribution [Eq. |(3J)] we then have v' c (t) = 0(v 2 \ A ), and 



from the approximate cluster expansion of Eq. ((5j> , 

Wcpmg(t) = cxp {-0(v 2 \ A )). 



(6) 



The analogous expression for the Hahn echo has 0(A 2 ); 
hence the extra symmetry of the time sequence in even- 
pulsed CPMG echoes leads to a cancellation of the lowest 
A 2 order. For this reason, we predict an enhancement of 
even CPMG echoes over odd echoes when this lowest 
order result is valid (we find that it is generally valid as 
long as the cluster expansion itself is convergent as we 
observed and explained in our Hahn echo analysis [6J]). 
Furthermore, Eq. (|6|) shows that the logarithm of the 
CPMG echo as a function of r (but not t) scales with the 
number of puses squared. 

Similarly, the lowest order with respect to the time 
expansion gives vc(t) = 1 — 0(t 6 ) so that i>cpmg( t ) = 
cxp (-0(t 6 )) for the short time behavior of the CPMG 
echoes. The analogous expression for the Hahn echo has 
0(t 4 ) which is cancelled [much like 0(A 2 ) of the inter- 
bath perturbation] by symmetry in the time sequence of 
even-pulse CPMG echoes. This short time behavior is 
exhibited by the GaAs system but not the Si:P system 
(the reason relates to the shape of their electron wave- 
functions and resulting cluster contribution statistics). 

Applying the lowest order inter-bath perturbation rep- 
resented by Eq. ([6]) to the GaAs quantum dot system 
of Fig. [1] and fitting to these numerical results yields 
ln(u CPMG (r)) » — v 2 (r/55 fis) —v 2 (r/31/zs) where the 
first (second) term result from clusters of size two (three) . 
Interestingly, 3-cluster contributions dominate as a con- 
sequence of the low order A 2 cancellation. This approxi- 
mation agrees with the 0(t 6 ) short time behavior that we 
anticipated. It also agrees with the exact results shown in 
Fig.fflexcept as shown by the dotted lines which deviate, 
in a conservative way (predicting overly quick decoher- 
ence), from their corresponding exact results. 

For the GaAs system of Fig. [IJ the two-pulse CPMG 
echo prediction of T 2 = 120 /zs [14j is more than four 
times the Hahn echo T 2 = 28 /xs, demonstrating the en- 
hancement of even echoes. More generally, in the con- 
servative perturbation approximation, v 2 r e factor im- 
plies that r effectively scales as v~ x l 3 \ this means that 
the time between pulses, 2t, must be shortened as one 
increases the number of pulses, 2v, in order to yield 
the same degree of coherence. However, since the to- 
tal pulse sequence time, t, is given by t = Avt, we have 
T 2 = v 2 ! 3 x 120 (is. Therefore, coherence is enhanced 
by applying multiple CPMG pulses, but one is experi- 
mentally limited by the minimum time allowed between 
pulses which is in turn limited by the time needed to 
apply each pulse. 

Applying the lowest order inter-bath perturbation rep- 
resented by Eq. ^ to the Si:P system of Fig. [3] (with 
B || [100]) and fitting to these numerical results yields 
ln( UcPMG (r)) » -„ 2 / 2 (r/230 ^s) 4 ' 3 - v 2 f {r/7A M s) 4 ' 3 



where / represents the fraction of Si that are the 29 Si 
isotope, the only isotope of Si with a net spin that can 
contribute to the spectral diffusion. The f 2 (/ 3 ) depen- 
dent term comes from 2-cluster (3-cluster) contributions. 
Note that isotopic purification will make the most impact 
when 3-cluster contributions dominate over 2-cluster con- 
tributions because of their relative scaling with /. These 
two contributions become comparable when / « 0.7%. 
This may be an important consideration for a cost-benefit 
analysis of isotopic purification of Si. 

For natural Si (/ = 4.67%), ln(u CPMG (r)) w 
—v 2 (r/0.63 ms) ' . This agrees with the exact results 
shown in Fig. [2] except as shown by the dotted lines 
which deviate, in a conservative way (predicting overly 
quick decoherence), from their corresponding exact re- 
sults. For the two-pulse CPMG echo (y = 1) of natural 
Si, T 2 — 2.5 ms, more than four times longer than the 
Hahn echo T 2 — 0.6 ms, demonstrating, again, the co- 
herence enhancement of even echoes over odd echoes. In 
general, the above i/ 2 r 43 factor implies that r effective 
scales as v~ 0A7 ; since t = 4vt, T 2 scales as v 053 . This is 
comparable to, but not quite as good as, the ;/ - 67 scal- 
ing of T 2 for GaAs. For the natural Si system of Fig. [2j 
T 2 = u 053 x 2.5 ms, in this conservative perturbation 
approximation. 

We have developed a cluster expansion theory for cal- 
culating the enhancement of spin coherence in the con- 
text of multiple-pulse spin echoes in the limit of large 
applied magnetic fields. We find considerable enhance- 
ment of spin coherence in semiconductor qubits when 
even CPMG pulses are employed; T 2 increases four fold in 
both GaAs and Si:P systems in comparing 1- (Hahn) and 
2-pulse sequences, and increases with additional pulses. 
Using a conservative approximation for these GaAs and 
Si:P systems, T 2 scales with the number of pulses to 
the power of 0.67 and 0.53 respectively. We also re- 
port a cross-over in the scaling of even CPMG echoes as 
a function of isotopic purification of Si in Si:P marking 
a transition between 3-cluster and 2-cluster dominance. 
This is a further consequence of the interesting symme- 
try related cancellation responsible for our important re- 
sult that even-numbered pulse sequences are subtantially 
more effective in extending/restoring spin coherence than 
odd sequences. 
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